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Abstract 


Assuming a stationary, radial, spherically symmetric solar wind and a 
radial magnetic field direction in the vicinity of the sun, an equation of 
motion for ions heavier than protons in the solar wind is derived. The ge- 
neral properties of this equation are discussed and the results of numeri- 
cal integrations are given. These results are based on the assumption of 
maxwellian velocity distribution functions for electrons, protons and ions, 
but the effects of first order deviations from such distributions are also 
presented and discussed. It is shown that dynamical friction, i. e. momen- 
tum transfer from protons to heavier ions accounts for the observed fact 
that heavier ions - if accelerated at all - normally reach the same velocity 
as the protons in the solar wind. Because of the non-linear relation bet- 
ween dynamical friction and proton-ion velocity difference a minimum pro- 
ton flux is required to carry a certain ion species in die solar wind. For- 
mulae comparing the minimum fluxes for different ions are given. It is 
shown that elements up to and beyond iron will be carried along in the 
solar wind as long as helium is carried along. Substantial isotopic fractio- 
nation is possible, in particular in the case of helium. The effects of ion 
motion and escape on abundances in the corona and in the outer convective 
zone of the sun are discussed. % 


1 . INTRODUCTION 
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It has been established for several years (SNYDER and NEUGEBAUER, 

1964; WOLFE et al. , 1966; HUNDHAUSEN et al. , 1967) that 4 He ++ ions 
are present in the solar wind and that their relative abundance is highly 
variable, values for He/H between 0. 01 and 0. 25 having been observed 
(HUNDHAUSEN et al. , 1967). The average He/H ratio in the solar wind 
is about 0. 04 - 0. 05 (HUNDHAUSEN et al. , 1967; SNYDER and NEUGE- 
BAUER, 1966; OGILVIE and WILKERSON, 1969; ROBBINS et al. , 1969), 
which is lower than the value of 0. 1 often assumed for the sun (cf. BISWAS 
and FICHTEL, 1964). At least during stationary conditions the bulk velo- 
cities of H and He are normally the same within a few percent (HUNDHAU- 
SEN et al. , 1967; OGILVIE et al. , 1968; ROBBINS et al. , 1969). Recent- 

3 

ly, also the presence of He and oxygen has been reported in the solar 
wind (BAME et al. , 1968). However, as these ion species are observed 
only under very favourable conditions, no figures on their abundance ave- 
rages and variations can be given. Abundances of still heavier elements 
such as Ca, Fe, and Ni have been determined in the corona (POTTASH, 

1963), and it appears that these elements are overabundant relative to the 
photosphere. 

Both the coronal and the solar wind observations indicate that there are 
processes at work which are substantially influencing the relative abundan- 
ces of ions. If these processes would be sufficiently understood abundance 
measurements in the solar wind could be used to gauge conditions in the 
solar atmosphere. A theoretical understanding of enrichment processes 
is also important if solar abundances, for instance of isotopes, are to be 
inferred from solar wind observations. 

PARKER (1961) discussed the static non- convective atmosphere formed by 
a multi- component plasma in a gravitational field. He showed that a mass 
separation of ions will occur similar but not identical to the mass separa- 
tion in an atmosphere composed of neutral molecules. Ions with mass/charge 
ratios larger than that of protons will accumulate at the bottom. This re- 
sult has sometimes led to the assumption that the abundance of elements of 
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higher mass number may be significantly reduced in the solar wind. 

JOKIPII (1965) and DELACHE (1965) used the diffusion equation as given 
by CHAPMAN and COWLING (1958) for studying the motion of ions in the 
corona. The interesting result is that in the lower corona heavier ions 
move markedly slower than protons, producing enhanced number densi- 
ties of the former. More recently DELACHE (1967) and NAKADA ( 1969) 
have extended these studies and have shown that essentially as a result 
of thermal diffusion, very strong enhancements of heavier elements may 
occur in the lower corona, unless convective mixing is dominant. Con- 
ditions for the occurence or absence of radial convection in the lower 
corona have been discussed by NAKADA (1969). 

The diffusion equation as given by CHAPMAN and COWLING (1958) is an 
approximation applicable only as long as bulk velocities and their diffe- 
rences are small compared to thermal velocities. Thus its applicability is 
restricted to the lower corona, since in the coronal expansion protons and 
in particular heavier ions reach high Mach numbers. 

For the purpose of theoretical study, two fundamental problems may be 
distinguished: 

1. Do there exist diffusion processes or perhaps other processes 
across the photo sphere -corona transition region which are en- 
hancing the abundance of heavier elements in the corona? 

2. Which are the factors determining acceleration and escape of 
ions other than hydrogen from the lower corona into interpla- 
netary space? 

In this paper we address ourselves mainly to the second problem. We 
shall derive and discuss a general equation of motion for ions applicable 
in the solar wind acceleration region. 
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2. THE PROTON" ELECTRON GAS 

The theory of expansion of the solar corona and of acceleration of the 
plasma to supersonic velocities was developed by PARKER (1958, 1962, 
1963, 1964). We shall summarize here those of Parker's results which 
are of direct relevance to the discussion in this paper. The equation of 
stationary radial motion is given by 

V sLL _L_ sL (/uk~) = c (2.1) 

ctr u r *i it ' / r* 

V is the bulk velocity of the plasma, n is the number density of protons 
which is equal to that of electrons. The conservation of mass can be 
written as 


«Vr‘ - (2.2) 

where is the proton flux at = 1 a.u. For radial flow s is equal 

to 2, whereas s > 2 corresponds to a flow diverging more rapidly. 

In his original work, PARKER (1958, 1963) used a polytropic law to com- 
plete die system of equations for the unknowns n(r), V(r) and T(r): 

T » Court. * (2. 3) 


By proper choice of a, heat transport can be taken into account. Later 
PARKER (1962, 1964), NOBLE and SCARF (1963) and WHANG and CHANG 
(1965) solved both the momentum equation and the energy equation by assu- 
ming that heat transport is solely conductive. More recently STURROCK 
and HARTLE (1966) pointed out that energy partition times between elec- 
trons and protons under coronal conditions are large enough as to give 
substantial temperature differences between these two components. These 
authors have accounted for this in their two-fluid model (cf. also HARTLE 
and STURROCK, 1968). For the purpose of this paper, we consider the 
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poiytropic model adequate, because the refinements mentioned can be 
approximated by adjusting a, or even by assuming it to be varying with 
solar distance. In fact the polytropic model is very convenient for dis- 
cussing the factors determining ion acceleration and therefore, our 
discussions will be based on it. 

The polytropic law and the continuity equation for the mass imply 
nkT o C (V* r 8 ) °. Inserting this relation in the differential equation (2. 1) 
we obtain 

where the temperature may be expressed in terms of the velocity as 
T oC (V* r*) 1 " 0 . Integrating (2.4) we arrive at the Bernoulli equation 

* 

is*!!* hi IcT _ „ f . (2.5) 

2 - r 

In order to determine the admissible values of the constant of integration 
E we observe with Parker that the differential equation (2.4) is singular 
if the velocity V reaches the critical value 

A. (M r V* sr* UT (2.6) 

As this point the gradient dV/dr becomes infinite unless the right hand 
side in (2. 4) vanishes simultaneously 

2 ntlcT * . (2.7) 

r 

The two relations (2.6) and (2. 7) define the so-called critical point 

(r , V , T ). Parker has shown that die physical solution must pass 
c c c 

through the critical point if the boundary conditions at r-+ o and at r-»o» 



are to be satisfied, i.e. 
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v&;- v t -(<3A/ e /r/ij ,/v 

7 ( 2 . 8 ) 

Tfoj - T« - <SH s U4 r /UMkr c 

This implies that only one particular value of the integration constant E 
is acceptable. Evaluating (2. 5) at the critical point one finds 

f • < 2 ’> 

To obtain a convenient representation of the physical solution we express 
it in terms of the dimensionless quantities: 

f - n. f 

y(r) - X. W(t) ( 2 . 10 ) 

T(r) •T,X(j) 

The dimensionless temperature T is given by 

, ( 2 . 11 ) 


and the velocity W is the monotonically rising solution of the transcenden- 
tal equation 




A 

•<-<f 


C . (2.12) 


The functions W( £ ) and Z ( J ) are determined uniquely once the poly- 
trope index a and die divergence coefficient s have been chosen. With 
W( J ) and T ( J ) known, the physical solution V(r), T(r) is given by (2. 10), 
(2. 8) in terms of die position of die critical point which is the essential 
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parameter in our representation of the physical solutions. 

For the particular case of isothermal expansion characterised by o -* 1 
(2. 11) and (2. 12) reduce to 


T~A 


*£ _ z*,(W f) 



( 2 . 13 ) 


The following asymptotic relations derived from (2. 12) and (2. 13) will 
be needed in our further discussions. 


*•> ± tf i . 

w-f^r* 4 f p + o(v] 

T - cfr-*) f* 

•«— * < 2 - ,6 > 

The condition ftat W should not go to infinity for ^ 0 places a limit 

on OS / m A 

(«<t for *m) u-17) 




(2- »*> 
(2. 15) 


b.) Pv>4 

w m Wo* -£5 (*£•?} * "f 


Wo. » l % (*' A ) 

is real only for 




•<< 




U-A 

The temperature decrease is given by 

r«c f 

For tea limiting case 0*1 


(•*< £ t mXj 




(2. 18) 

(2. 19) 

( 2 . 20 ) 
( 2 . 21 ) 
( 2 . 22 ) 


is obtained. 



3. EQUATION OF MOTION FOR IONS 


Before we derive an equation of motion for ions in the solar wind a dis- 
cussion of the role of magnetic fields is called for. The question is, 
whether heavier ions can be swept along by the magnetic field lines 
moving with the proton-electron gas. It should be noted that in order 
to move with the bulk velocity in the solar wind, heavier ions have to 
acquire a kinetic energy many times their original thermal energy. It 
is often assumed that the lines of force near the sun in the acceleration 
region of the solar wind are approximately radial (cf. PARKER, 1963; 
JOKIPII, 1965). This assumption is probably justified under quiet-day 
conditions in most of this region. In fact PARKER's (1963) "spiral mo- 
del" of the interplanetary magnetic field predicts nearly radial field 
lines in the vicinity of the sun. A radial field, even when divergent, 
would of course not give an effective radial acceleration. In an exactly 
transverse field moving with the plasma any ion would have the same mean 
radial velocity as the plasma, irrespective of its initial velocity. However, 
in the case of oblique field lines, neglecting collisions the mean radial 
velocity of an ion would depend on the initial velocity. If for instance an 
ion is placed into the moving plasma with an initial velocity around aero 
the radial motion of the ion is slower than that of the plasma and non- ra- 
dial ion motion occurs. Thus an ordered magnetic field of spiral configu- 
ration cannot accelerate ions to full solar wind velocity. This, however, 
would be in contradiction with the observation that H and He have virtually 
identical mean radial velocities (HUNDHAUSEN et al. , 1967; OGILVIE 
et al. , 1968; ROBBINS et al. , 1969). A highly disordered magnetic field 
would be more efficient in acceleration, but still the magnetic field as the 
prime cause of acceleration appears to be insufficient, because the He and 
H velocities are so permanently and closely equal. The conclusion is, that 
much of the salient features of ion motion in ‘’he quiet-day solar wind should 
be found in a model neglecting the effects of the magnetic field. 

We assume a plasma essentially composed of electrons and protons with a 
small abundance of another ion species in addition. The state of the plasma 
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is described by the partition functions of electrons, protons and ions, 

f (5f, t), f (5f, T, t) and f_ A (1T, t) respectively. The evolution in 

e p Zt, A 

time of these partition functions is described by the Boltzmann equa- 
tion whose explicit form in the present situation is given in the appendix. 
To simplify the problem, we assume that each of the components of the 
plasma is in an approximate local thermodynamic equilibrium, i.e. , 
we assume that the velocities of the various components V = (e; p; Z, A) 
are distributed according to Maxwell's law around the individual bulk 
velocities v^ (ST, t) with a spread given by the temperature T ^ (x, t): 

- •vf-SF’-V/*** / (3 ‘ 11 

The quantity n ^ (3fr, t) denotes the particle density of the component V . 

In this approximation thermal diffusion is neglected. Thermal diffusion 
is proportional to the temperature gradient and - as pointed out by 
JOKIPII (1965) - its effects are relatively small in the corona outside 
the sharp temperature rise in the very narrow photosphere corona tran- 
sition region which is not included in the quantitative theoretical treatment 
of this paper. However, we shall see in Chapter 9 that also in the corona 
thermal diffusion may have an effect on the motion of ions. 

In order that representation (3. 1) of the partition functions be an approxi- 
mate solution of the Boltzmann equation, the quantities n v , v v , and T t , 
which describe the state of the plasma in this approximation, must satis- 
fy the transport equations for particle number, momentum and energy 
respectively. The transport equation for the particle number is the fami- 
liar conservation law 

( 3 - 2 ) 

The analogous differential equation describing transport of momentum, 
i.e. , the differential equation for "v ^ (x“, t) is derived in the appendix. 
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The transport equation for Ty (x, t) is of a similar nature. In the pre- 
sent paper we do not attempt to solve this system of differential equa- 
tions for n y , ^ , and Ty . Instead, we simplify the problem further 
by assuming that the energy transfer between the different plasma com- 
ponents is sufficiently effective such that to a good approximation the 
plasma is characterised by a single local temperature T^T^T^ ^*T. 
Actually, as mentioned in Chapter 2, a significant difference between 
electron and proton temperatures in the solar wind cannot be ruled out 
and the Parker model has been generalised by STURROCK and HARTLE 
(1966), to account for this difference. There is no reason, however, to 
expect that this modification strongly affects the motion of the ions. 
Concerning the equality of proton and ion temperatures we refer the rea- 
der to SPITZER's (1962) relation for equipartition times between protons 
and ions and conclude that under coronal conditions the ion temperature 
is indeed closely bound to the proton temperature. 

To summarise this discussion, we list the approximations involved in 
our treatment of ion transport: 

a) Forces die to the magnetic field are neglected. 

b) The velocity distribution of each particle species Electrons, 
protons, ions) is maxwellian around an individual bulk velocity. 

c) The local temperatures of electrons, protons and ions are equal. 

d) nj <c n p * n # « n 

In this approximation the equations for stationary motion of one particu- 
lar ion species are 

*1 >fr m CCiu*. 

where v denotes the bulk velocity, n^ the density, Z e the charge and 

Am the mass of (he *on species under consideration. The explicit 
P 

expression for die collision term C is (cf. BURGERS, I960). 

(3.4) 


U ItT A •) 
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A derivation is given in the appendix. In corona and interplanetary plasma 
JLa A i* nearly constant, A* A te 22 (PARKER, 1963). 

The function C (x) was first introduced by CHANDRASEKHAR (1943) and 
is tabulated by SPITZER (1962). It is defined as 

2 * k (3 - 5) 

4>(x) being the error function. The asymptotic behaviour 

G (x) = 0. 37 x for x« 1 
G (x) = l/Z for x fcV 1 

shows that "dynamical friction" (CHANDRASEKHAR, 1943), i.e. accele- 
ration due to collisions is proportional to the difference in bulk veloci- 
ties until V-v approaches the thermal velocity of the protons, beyond 
which it is actually decreasing. 

Equation (3.4) includes only collisions by protons. The momentum trans- 
fer by electrons is smaller by a factor 

(m p /m e )^ = 43 for |V-v| « (2KT/m p )^ and T e fc T . The effects of 
electrons become dominant only when die argument of G in (3.4) is larger 
than 3. 8. 

If we neglect the perturbation of the proton and electron motion due to 

collisions with ions and use n = n = n, v = v = V then the transport 

e p e p 

equations for protons and electrons take die form 



±*-(MkT)+<.E «o 


In the electron transport equation we have neglected the kinetic term 

s'-VdV/dr since it is of order m /m as compared to the remainder. 

e p 

Elimination of the electrical field E in (3. 6) leads to the basic plasma 
equation (2. 1). On the other hand, we may use (3. 6) to determine the 
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electric field 



±- f §-& + 

a- l r v 



(3.7) 


Inserting this expression in the ion equation (3. 3) we obtain 



This equation, supplied with suitable boundary conditions at r — and 
at r — , determines the bulk velocity of the ions, once the proton 
velocity V(r) and the temperature T(r) are known from a model of the 
main solar wind components. In the following we restrict ourselves to 
Parker's polytropic model for the proton- electron plasma. We use the 
representation (2. 10) of this model and introduce the dimensionless ion 
velocity w(^ ) as 


y(f) - V t w (fj ( 3 . 9) 

where is file velocity of file proton gas at the critical point. In terms 
of this variable file differential equation (3. 8) may be rewritten as 


^f( 2w -Jk 

The dimensionless quantity f denotes the reduced flux 


(3. 10) 




?k-t 


(3.11) 
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where is a convenient flux unit 

4> (.HI t* A.) *uy rl /l ( GMq/i f l 

m A-3 « A<f Om x it c'" 1 


(3. 12) 


and the symbol 4W stands for essentially die ratio of the velocity 
difference to the thermal velocity 



Note that the reduced flux f depends on the critical radius of the proton 
gas and thus on the temperature in the corona. 

Equation (3. 10) represents the basic equation for ion transport in the 
framework of Parker's polytropic model of the proton- electron plasma; 
the following chapters contain a discussion of the properties of this 
equation. 
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4. DISCUSSION OF EQUATION 3. 10 FOR LIMITING CASE OF f=0 

In this Chapter we shall discuss the limit £ — ► 0, in which case the 
collision term in equation (3. 10) becomes negligible. (3. 10) is re- 
duced to 


if %(\) (4 - i) 

with 

W( | ) and T (f ) are known from the polytropic model, thus g is a 
known function of J , and (4. 1) can be integrated. w(^ ) determined in 
this way includes one free integration constant. 

a) Isothermal expansion (a = 1). 

As an illustration we shall first discuss the isothermal case. Here 
equation (4. 1) can be integrated in closed form: 

w*-- jL. w ^ i±A w*- w + JL. (“-3) 

A ZA A A $ A 

K is the integration constant. Generally for a given K, there exist two 
values of w for each ^ , thus the solution is split into an upper branch 
w + and a lower branch w . Physical boundary conditions are determi- 
ning the choice of K: w should be limited for £ « 1, and n^ should go to 
aero for | . In order to apply these conditions, die asymptotic be- 

haviour of w( j ) has to be analyzed. From (4. 3) we obtain for the upper 
and lower branches for f « 1 
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W 4 ©C 


f 




(4.4) 


w. = f - K+r-d j 


of which only the lower branch is acceptable. 
For ^ —»oo we obtain 






(4. 5) 


(4.6) 


»- « ("“"Vi. 




(4.7) 


From n^-e 0 follows (w jf 1 ' )— *oo , a condition which is fulfilled only by 
the upper branch. Thus the solution of physical interest has to cross 
over from the lower branch at ^ « 1 to the upper at ^ -»** and there 
exists only one solution, the critical solution, with one specific value 
of K, passing through the critical point. At this point (^ c t W c ) the left 
and right sides of (4. 1) vanish simultaneously, yielding the relations 






.**/» 


(4.8) 


0 

t _ Jiii dzX-di i 
> c 




( 4 . 9 ) 


(4. 10) 


The critical point for protons falls of course on W( ^ ), but it is not iden- 
tical with the critical point in Parker's theory of the proton- electron gas 
(W * 1, £ = 1). There is no contradiction in this result, because the 

only solution of physical significance of the families of curves defined by 
(4. 3) is the critical solution, and for protons this is identical with W(^ ) 
as given in (2. 13). 


* 
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We conclude that for isothermal expansion, heavier ions are accelerated 
to supersonic velocities. In the case of low solar wind flux their veloci- 
ties remain, however, lower than the solar wind velocity W. From (2.22) 
and (4. 6) we obtain 


w 


(*af 


Formally there is no limit for A. However, for A< 


(4.11) 

and Z/A — O 


the critical point 


h 


goes to infinity. 




14. 

i 


(4. 12) 


i. e. , it moves to a solar distance where an isothermal approximation 
becomes meaningless. 

b) Polytropic expansion (a >1). 

The question to be discussed is whether equation (4. 1) has a solution for 

which w is monotonously rising with j: . In this case, the left side of 

(4. 1) goes through zero for a critical velocity w which is reached at 

c 

critical point ^ ^ : 

<413) 


g(f ) as defined by (4. 2) vanishes at ^ . In order to establish the 
existence of a monotonously rising solution w ( ^ ) , it is sufficient to show 
that g(^ ) is negative for ^ c , and positive for p > ^ , and that 
solutions w(^ ) exist which meet the physical boundary conditions. As in 
the isothermal case we have generally two branches w + ( | ) and w_(^ ) 
for a single integration constant. From (4. 2) and the approximations 
given in Chapter 2 we obtain for £ <* 1 : 


J. .2 A-i-(t+4)*ad 


( 4 . 14 ) 
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For (2 A-Z)»l, g Q it positive and g - as required - negative, as long 
as 


w< (xmJ /r 


(^<1 for t-Z) (4. 15) 


The upper and lower branches of w for t « 1 are: 


W t oC f 

W.*c 


(4. 16) 


The lower branch meets the boundary condition, because g Q > 0. 
For £ » 1 equation (4. 2) gives 


ffi) •*(**'•) f 



(4. 17) 


g( ^ ) is positive, again as long as condition (3. 15) for a is fulfilled. 
The upper and lower branches of w( ^ ) for ^ » 1 are: 




— IVos ^ ' 


Wqo 2A 


»_-C f 




(4. 18) 


Like in the isothermal case the lower branch is not compatible with 

n. -♦ 0, and the upper branch is the physically acceptable solution, as 
1 / 1 \ 

^ " a " . We conclude that a monotonously rising solution w( J ) 
exists as long as a< (s+l)fc. This condition is identical with the limit 
placed on a for polytropic expansion of the proton- electron gas (cf . 2. 17). 


5. HIGH SOLAR WIND FLUX 


If we let the tolar wind flux increase, dynamical friction, i. e. the colli- 
sion term in equation (3. 10) forces w closer to W. In the limiting case of 
very high fluxes (f »1), W-w will be small compared to W, and w can be 
approximated by W in equation (3. 10) everywhere but in the collision term: 


4 W • ^ 



(5.1) 


where the function 



is known from the proton- electron model: 





(5.2) 


In this approximation the original differential equation for the ion velo- 
city degenerates into an algebraic equation for which no boundary condi- 
tions have to be supplied. Because equation (5. 1) is based on the assump- 
tion of maxwellian velocity distributions it does not contain a thermal 
diffusion term. In this respect it deviates from the diffusion equation as 
given by CHAPMAN and COWLING (1958). The consequences of this 
deviation will be discussed in Chapter 9. 
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6. DISCUSSION OF THE GENERAL SOLUTION 
The equation of motion (3. 10) can be written as: 


( 6 . 1 ) 

^ ZW ~ Jt* 

g(^ ) is defined by (4. 2), and AW by (3. 13). The singularities of 
equation (6. 1) (simultaneously vanishing nominator and denominator on 
the right side), determine the topology of the family of solutions. In 
some cases, one such singularity exists, and the corresponding family 
of solutions is sketched in Figure 1. In other cases there are three 
singular points, and the family of solutions has a topology as shown in 
either Figure 2 or Figure 3. In any case only the one solution passing 
through the critical point P is compatible with the physical boundary 

G 

conditions. 


The critical velocity w is given by equation (4. 13). Inserting w = w 

G C 

in the collision term, can be determined. If f is increased P^ moves 
towards the left along the curve w = |T( j: )/2 w A| ^ . Thus the colli- 
sion term has the general effect of lowering and bringing w closer 
to W. 


However, even for strong fluxes the velocity of the ions does not necessa- 
rily approach die bulk velocity of the solar wind for large solar distances. 
In the limit ^ 0* there are two conditions to be met, if w is to approach 

W. 


1 . ) For the collision term C must dominate die function g( ^ ) . 

Using die approximations (2. 19), (2. 21), and the fact that die function G 
is limited we obtain 


c < 


f 




( 6 . 2 ) 
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Comparing this with the asymptotic behaviour of g(^ ) as given in (4. 17), 
we obtain 

> l***£d. £ld f j jttfrxj (6.3) 

2.) In order to have an effective collision term, C must rise with rising 
W-w (cf. 3i9 s s 4 3. 5). This is assured only if W-w approaches aero 

Vz 

faster than X , in order that AW (cf. 3. 13) remains smaller than 
unity. 

We can study the consequences of this condition by solving the equation 
of motion (3. 10) for large solar distances. If we assume an asymptotic 

Vz 

behaviour W-w < T , and introduce 

S * W-w (4.4) 


we obtain (cf. 5. 1, 5. 2 a n d 6. 3 ) 


For f » 1 and 4W «C i/X 1 (6. 5) is approximated by (cf. 2. 21) 



a and b are positive constants, and 


( 6 . 6 ) 


M 



-J 




Ut 


(6.7) 
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It it easily shown that this first order linear equation has a solution 
with 


.&U4 J » O 

only when n< 0. From (6. 7) follows 


( 6 .*) 


cx: > A + A- Mr t.l) (6.9) 

It can be shown that condition (6. 9) is applicable not only for PARKER's 

(1963) solar wind model but also for the two-fluid model of STURROCK 

and HARTLE (1966), if a is replaced by a , the effective polytropic 

P 

index of die protons. According to (6. 9) the natural increase of a (res- 
pectively a ) with larger solar distances (cf. PARKER, 1963; STURROCK 
P 

and HARTLE, 1966), will force heavier ions to travel with the velocity 
W of the proton- electron gas, unless somewhere at intermediate distances 
the velocity deficit W-w was large enough to make AW 1, in which case 
the collision term became ineffective (cf. 3. 5). 


7. ESTIMATE OF "MINIMUM FLUXES'* 


In this Chapter we shall discuss the general conditions under which heavier 
ions will travel with the bulk velocity in the solar wind. It is assumed that 
a increases gradually from a value 'ose to 1 near the sun until it approaches 
a value around 1.4 to 1.5 (cf. STURROCK and HARTLE, 1966). At a cer- 
tain solar distance jp the value a = 4/3 (6.9) will be surpassed. If, because 
of insufficient flux strength, W-w has become so large at ^ that dw > 1, 
then the collision term becomes ineffective beyond and the ion velocity 

w will not converge towards the bulk velocity W. A minimum flux f^ 
sufficient to bring the velocity of a given ion species up to the solar wind 
velocity can be calculated by numerical integration of equation (3. 10). The 

actual values of f . will of course depend on the details of the chosen solar 
min 

wind model, and in particular the choice of s will greatly affect the minimum 
flux f . However, the relative magnitudes of different ion spe- 

cies will depend much less on such details. 

As shown above an approximate condition for obtaining ions traveling with 
W is 


If we use the high flux approximation (5.1), we obtain 


(7.1) 



(7.2) 


Because F ( ^ ^) is model dependent, equation (7. 2) does not readily yield 
absolute values fcr minimum fluxes. However, (7.2) should be a valid 
approximation if used for calculating ratios of minimum fluxes of different 
ion species: 
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In the lower corona the times for an atom to approach equilibrium ioni- 
zation are short compared to the characteristic times of bulk motion 
(JOKIPII, 1965). Hence it can be assumed that all ions have approxima- 
tely attained an equilibrium ionization corresponding to the temperature 

in the lower corona. For a number of elements and isotopes, these de- 

6 o 

grees of ionization are given in Table 1 assuming = 10 K. They were 

calculated by the method of ELWERT (1952) (cf. BILLINGS, 1966) using 

the ionization potentials given by LOTZ (1967). Also given is the factor 

A) which in the approximation discussed here determines the 

minimum flux. According to Table 1 isotopic discrimination in solar wind 

3 , 

acceleration can be appreciable for light elements. Particularly the He /He 

3 (l) 

ratio could be changed in favour of He , because the P' ' values of these 

two isotopes differ by a factor of 1.7. The important and perhaps somewhat 

surprising conclusion is that if medium and heavy elements are present 

in the lower corona, they are accelerated to solar wind velocity, as long 
4 

as He is accelerated. With the exception of deuterium, which should not 

4 

exist in the sun anyway. He is no more likely to be accelerated to solar 
wind velocity than all the stable nuclear species up to and beyond iron. We 
have not exactly calculated P^ for heavy elements, because LOTZ's (1961 
table of ionization potentials ends with zinc. However, estimates based 
on the Thomas-Fermi model indicate that P*^ for very heavy elements 
will not be higher than 1. 25, the value for helium. 

If, as asserted in this Chapter, a minimum proton flux is necessary to 
carry an ion species in the solar wind, one would expect to observe a corre- 
lation between the He/H ratio in the solar wind and the flux. Such a corre- 
lation has not been found (ROBBINS et al. , 1969) in the Los Alamos Vela 3 
data (cf. HUNDHAUSEN et al. , 1967). It should be remembered, however, 
that the minimum flux is dependent on a number of parameters such as 
temperature in the corona, effective radial expansion index s, or polytropic 
index oc . Moreover, the model presented here is stationary, and non- 
stationary processes are likely to affect the He/H ratio. The He abundance 
in the solar wind may in part also reflect the variable influence of thermal 
diffusion between photosphere and coronal base. All these effects would 
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contribute in masking a correlation between abundance of one ion species 
and proton flux. However, these effects should cancel out to a large extend 
when the behaviour of different ion species is compared. 
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8. RESULTS OF NUMERICAL INTEGRATION 


Numerical integration! of equation (3. 10) were performed for • s 2 
and uniform polytropic index a, choosing a = 1. 1 and 1.4. Actually, for 
quiet solar conditions the effective value of a will vary from about 1. 1 
near the sun to about 1.4 to 1.5 at larger solar distances (cf. the 
values for protons in the two-fluid model of STURROCK and HARTLE, 
1966, 1968). Thus, our results with a = 1. 1 should give a realistic pic- 
ture of the physical situation near die sun, while the solutions with a * 1.4 
should be a better approximation for larger solar distances. 


4 2+ 20 8+ 56 14+ 

In Figure 4 the flux values for He , Ne and Fe and a = 1.4 are 
plotted versus the corresponding critical radii ^ The significance of this 
plot will be explained using neon as an example. For f = 0 the solutions 
have the topology of Figure 1, and ^ ^ is very high (3450). With increasing 
f the critical point moves continuously towards the sun until the flux f reaches 
a value of 3.4. For 3.4 < f < 4. 8 three singular points exist for a given f. 

If f is increased above 3. 4 the topology of the solutions is initially of the 
type shown in Figure 2, i. e. the critical point remains at very high solar 

distances £ . Somewhere between f = 3.4 and f = 4. 8 the topology changes 

| *»» . * 

to that given in Figure 3. Plots of the direction field of the differential 
equation and the numerical integrations show that this always happens for 
a flux which is only slightly above the relative minimum in the f - ^ ^ plot. 
Thus die critical point ^ ^ jumps discontinuously from a very high value to 
a place slightly to the left of the relative minimum in the f - ^ plot, in the 

case of neon from 200 down to ^ d 2. Figure 4 clearly demonstrates 

the idea of a minimum flux for a given ion species. Proton fluxes leading 
to a ^ c on the right side of the discontinuity are not able to carry the ions 
in the solar wind stream. If the flux is high enough to place the critical 
point to the left of the discontinuity (left of the relative minimum in die f - p 
plot) the collision term in equation (3. 10) becomes effective and pulls the 
ion velocity w towards the proton velocity W. The flux value at the relative 
minimum in the f - ^ c plot is therefore a good approximation for the mini- 
mum flux f ^ n defined in Chapter 7. This is born out by the numerical inte- 
grations. Thus for a: 1.4 minimum fluxes can be read off in good approxi- 
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mation directly from the relative minimum in the f - is plot (Figure 4). 

V c 

When the proton flux is below this minimum value, the corresponding ions 
are virtually not accelerated. If the flux reaches the minimum value the 
ions are suddenly "locked* in" with the protons. 

The behaviour of the ions described here is related to the runaway pheno- 
menon of electrons in a strong electrical field (DREICER. 1958). As dis- 
cussed by TUCK (I960), ALLIS (I960) and BURGERS (I960) runaway of 
electrons occurs when dynamical friction decreases with increasing drift 
velocity. 

4 2+ 56 144- 

Figure 5 shows results of the numerical integration for He and Fe 
It is seen that the minimum flux for helium in this case is about 6. 7 where- 
as for iron it is about 3. 8. Thus, for the a - 1.4 expansion model the mini- 
mum flux for helium is 1 . 8 times higher than that for iron. Approximation 
(7. 3) (cf. Table 1) gives 2. 5 for the ratio of minimum fluxes for these two 
ion species^ in fair enough agreement with the numerical result. It should 
be reemphasised that we do not attach great significance to the absolute va- 
lues for minimum fluxes given here, because they are model dependent. 
Nevertheless, it is rewarding that physically reasonable numbers are ob- 
tained (cf. Caption of figure 5). 

Figure 6 gives the ion velocity w as a function of solar distance ^ for a 
few flux values with a = 1.1. Here, of course, critical fluxes cannot be 
established, because a does not rise above the value 1. 33 which is required 
if ion velocities are to converge towards the proton- electron velocity (6.9). 
Still it is seen that for a given flux the velocity of iron is much closer than 
the helium velocity to the proton-electron velocity. Roughly identical w( ^ ) 
are obtained if the flux is 2. 5 times higher in the case of helium than in the 
case of iron. This is in good agreement with equation (5. l). 

As pointed out above, a: 1.1 should be a good approximation for the lower 
corona. For T^ = 10^ °C, the solar rim is at ^ = 0.2. It is seen that near 
the sun ion velocities can be substantially lower than proton velocities. For 
instance for f = 5 the iron velocity near the sun is more than 10 times lower 
than tike proton velocity. If a were to rise from 1. 1 to about 1.4 in the range 
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0.4^ £ 0. 8, at this flux the iron velocity would still converge towards 

the solar wind bulk velocity. In this case we would have an enhancement 
of a factor of 5 to 10 in the number density of iron in the corona as com- 
pared to the iron flux abundance in the solar wind at 1 a.u. 
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9. INFLUENCE OF NON- MAXWELLIAN DISTRIBUTION FUNCTIONS 

The equation of motion for ion* (3. 3) is derived from the Boltzmann equa- 
tion under the assumption that the velocity distribution functions are lo- 
cally maxwellian, which is of course an approximation to the true distri- 
bution function. Better approximations of f(x, t) have to be calculated 
and then used in the pressure and collision terms of (3. 3) in order to 
arrive at an equation of motion of higher precision. If the mean free path 
X is short compared to the characteristic dimension L of the system (L is 
the length over which an appreciable change of say the temperature T(r) 
occurs, L T /~ ) , then the plasma is in an approximate local equilibrium, 
i. e. the distribution functions of the plasma components approximately 
satisfy Maxwell's law with one and the same temperature and a unique bulk 
velocity. Of course in this approximation the dynamical friction vanishes. 
Short X implies a strong flux, and our high flux expression (5.4) for the 
difference of the bulk velocities of protons and ions represents a statement 
about the c&iation from local equilibrium. The order of magnitude of the 
velocity difference given by (5. 1) is V^'X/L as to be expected. This small 
effect is not necessarily exactly described by our simple method, which is 
based on the assumption that the velocity distribution of the individual 
plasma components satisfy Maxwell's law. In order to determine the velo- 
city difference correctly to first order in x/L we have to account for the 
corrections to the velocity distribution which arise in first order. These 
modifications may be worked out by expanding the Boltzmann equation in 
powers of x/L according to the Enskog- Chapman method given in detail by 
CHAPMAN and COWLING (1958). Using the general equation of these 
authors (8.4. 7) we arrive at the following expression for the first order 
difference between the bulk velocities of protons and ions ("diffusion equa- 
tion"): 


4W 


_d_ jlA-i-A . 


(9.1) 
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The thermal diffusion factor a^, (CHAPMAN, 1958) can be approximated 
by 

(9 ' 2) 

It is seen that our high flux approximation (5. 1) is equivalent to the Chap- 
man-Cowling diffusion equation, except for the thermal diffusion term. 

This correction results from the deviations from maxwellian velocity dis- 
tributions which were neglected in the derivation of (3. 10) and (5.1). 

In Chapter 7 we have used the high flux approximation (5. 1) in order to 

arrive at an expression for the minimum flux f (7. 3). Using the improved 

min 

high flux approximation (9. 1) in the same way we obtain 

M «’•*» 

In as much as (9. 1) is a better approximation than (5. 1) to the first order 
difference between the proton and ion velocities, equation (9. 3) should 
better than equation (7. 3) approximate the relation between minimum fluxes 
for different ion species. 

( 2 ) 

P '(Z, A) factors for several types of ions are listed in Table 1. o = 1. 25 
is chosen as a reasonable average. For Z = 1 the exact expression of a^, 
given by CHAPMAN (1958) has been used. It is seen that theft's vary 
less than the H 1 *'.. i... differences between minimum fluxes for different 
ions are somewhat smaller when using (9.4) instead of (7. 3). Nevertheless 
the general conclusions drawn in the proceeding chapters remain valid. 

Thus He* still has the highest minimum flux among all ions up to and be- 
yond krypton, except again for the heavy isotopes of hydrogen. Isotopes 
have appreciable differences in minimum fluxes. Thus isotopic fractiona- 
tion has to be expected. In particular die differences between die minimum 
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3 4 

fluxes of He and He are quite large and therefore considerable varia- 
3 4 

tions of the He /He ratio might occur in the solar wind. 

Clearly for those regions where the solar wind flux is strong in the sense 
that the mean free path X is indeed small compared to the characteristic 
geometrical length L the diffusion equation is the appropriate tool to in- 
vestigate the behaviour of the ion velocity. An analysis of this equation in 
the context of ion transport in the lower corona has been given by JOKIPII, 
(1965), DELACHE (1965, 1967) and NAKADA (1969). Our approach, on the 
other hand, provides a method to analyze the general case of an arbitrary 
solar wind flux. The basic assumption of our method, viz. the assumption 
that the velocity distribution of the individual plasma components may be 
approximated by a maxwellian distribution must be regarded as a pheno- 
menological parametrization of the real velocity distribution in terms of 
its center (bulk velocity) its height (particle density) and its width (tempe- 
rature) . Accordingly we expect oiir method to describe the main features 
of ion transport only; to improve the validity of our approximation we 
would have to use a more refined parametrization of the velocity distribu- 
tion which allows for deviations from Maxwell's law. 


- 30 - 


10. SUMMARY OF CONCLUSIONS 

The model discussed in the proceeding chapters which is restricted to 
stationary conditions and assumes spherically symmetric flow and ra- 
dial magnetic field lines near the sun, leads to the following conclusions: 

For a species of ions heavier than protons there exists a minimum solar 
wind flux above which the ion velocity is virtually equal to the proton- elec- 
tron velocity in the solar wind at large solar distances. A condition for 
this is, however, that a , the effective polytropic index for protons rises 
above 1. 33 with increasing solar distance. During quiet solar conditions 
Op will rise beyond this value at some distance from the sun (one to a few 
solar radii) because of slow heat transfer between electrons and protons 
(STURROCK and HARTLE, 1966). For stationary proton fluxes below the 
minimum value ions will normally form a static atmosphere. 

The minimum flux for an ion species depends on the temperature in the 
corona, the radial expansion index s, and the polytropic index a. 

If at times should remain below 1. 33 throughout the region effective in 
solar wind ion acceleration, ion velocities would not reach the proton velo- 
city. In this case, the effects of the interplanetary magnetic field moving 
with the proton- electron gas have to be considered. If this field has an 
ordered spiral configuration (PARKER, 1963), then ions travelling along 
the field lines with radial velocities smaller than the proton velocity might 
occur in interplanetary space. 

The critical fluxes for medium heavy elements are smaller than for helium. 
Thus under stationary conditions these elements should be present in the 
solar wind as long as helium is present. 

In file lower corona ions heavier than protons may travel considerably 
slower than protons, giving rise to enhanced number densities for these 
ions, compared to their abundance in the solar wind. In addition there 
exists the mechanism of an enrichment of heavier ions in the corona by 
thermal diffusion between photosphere and corona. This mechanism, if 
not obliterated by convection, would not only further enhance number 
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densities of heavier ions in the corona, but also increase the abundances 
of these ions in the solar wind flux escaping from the sun. It is important 
to distinguish these two possible mechanisms for the enrichment of hea- 
vier elements in the corona, because only the latter mechanism will de- 
plete heavier elements in the sun. Comparison of abundances in photos- 
phere and corona with those in the solar wind should bring about this 
distinction and decide, whether the solar wind could have significantly 
changed die abundances in the outer convective aone of the sun. 
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TABLE 1 


Minimum proton flux factors P for various ion species. Degree of 
ionisation Z is calculated after ELWERT (1952) and LOTZ (1967) 

(cf. BILLINGS, 1966); for Kr it is estimated. P^ is obtained from 
equation (7. 3), and P ' ' from (9.4). It is seen that in both appro- 
ximations minimum fluxes for heavier ions are not higher than for 


Nuclide 


Degree of 
Ionisation 
Z 


Minimum Flux Factor 

( 1 ) ( 2 ) 

P P (a=l. 25) 


H 1 

1 

H 2 

1 

H 3 

1 

He 3 

2 

He* 

2 

c 12 

5 

o 16 

6 

N. 20 

8 

22 

Ne 

8 

28 

Si 

10 

» 36 
A r 

9 

_ 56 
Fe 

14 

* 

00 

u 

14 


0.00 

0.00 

2. 00 

2. 19 

4.00 

4.28 

0.75 

1.06 

1.25 

1.61 

0.72 

1. 19 

0.69 

1. 17 

0.48 

0.97 

0.55 

1.05 

0.45 

0.95 

0.77 

1.28 

0.49 

1.01 

0.78 

1.30 


Figure 1. Topology of solutions of equation of motion for ions (3. 10) 
in the case of one singular point. ^ , W and w are defined by (2. 10) 
and (3.9). The heavy line crossing the critical point (singular point) 
is the only solution for w( ^ ) compatible with the boundary con- 
ditions . 

Figure 2. Topology of solutions of (3. 10) in the case of three singu- 
lar points. Here the critical point P is the singular point on the right, 

c 

and the only solution w( ^ ) compatible "dth the boundary conditions 
passes through P . 

Figure 3. Again three singular points, but now the critical point is 
on the left. The solutions change from the topology in Figure 2 to that 
in Figure 3 through increase of proton flux. 

Figure 4. Dependency of the position of the critical point p c on t * ie 
proton flux f (cf. 3. 11) for polytropic expansion with a = 1.4. As shown 
in Chapter 8, the flux at the relative minimum in the f - p c plot 
corresponds approximately to the minimum flux defined in Chapter 7. 

Figure 5. Proton velocity W( ) and ion velocities w( ^ ) for diffe- 
rent proton fluxes f and as 1.4. Proton fluxes in protons/cm^ sec at 
1 a.u. ( ^ = 51.5) are given by = 0. 34 x 10® f cm * sec * for 

T = 10® °K at ^ = 1. W, w = 1 corresponds to Yj V = 152 km/sec. 

Figure 6. Proton velocity W( f ) and ion velocities w( ^ ) for diffe- 

' *2 

rent proton fluxes f and a = 1.1. Proton fluxes in protons/cm sec at 

8 “2 • i 

1 a. u. ( ft = 40.5) are given by = 0.3 x 10 f cm sec for 

T = 10® °K at 6=1. W, w = 1 corresponds to Vj V = 135 km/sec. 
Under the same conditions the solar rim is at 6 = 0. 19. 
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Appendix 


B OLTZMANN EQUATION : SMALL ANGLE SCATTERING OF CHARGED 

PARTICLES 


In this Appendix we briefly sketch the derivation of the transport equation 
for the bulk velocity used in our discussion of ion transport. We consider 
a plasma consisting of a number of different components a. The state of 
the system is described by a set of distribution functions f^fxt if, t) 
satisfying the Boltsmann equation (cf. HUANG, 1963). 


* u ■ t ♦ j— *■< ’ j f* (*> M i*) m 5Z 


(A. 1) 


The quantity C .(at, if, t) represents the collision term for encounters bet* 
op 

ween particles of type o and 6 and is explicitly given by 


Cy(W)- JVv/r-7/^, Jttv; 

-i LFXvj, (*/,*){ 


(A. 2) 


where dff*^ denotes the differential cross section of the elastic process 
a (i!) + 0(v) — * o(u) + 3(v*) . The dominating contribution to the differential 
cross section arises from small angle coulomb scattering, described fay the 
Rutherford formula (c. m. system). 


^ (Vy * {k *-ir) dsi 


(A. J) 


In this expression 4 ^ stands for the charge of die particles of type a, 

n _ denotes the reduced mass 
o3 


/iu, - K •*/• / ♦ *|t 


(A. 4) 



and w it the relative velocity of the two particles, w = vf- v? If the Ruther- 
ford cross section is inserted in (A. 2) one obtains an integral that diverges 
at 0 » O . As is well-known the singularity of the Rutherford cross 
section at 0*0 reflects the long range nature of the coulomb force va- 
lid between isolated particles. In a plasma the long range tail is cut off by 
shielding effects. One accordingly cuts the integral over the center of 
mass scattering angle & off at a value ^ min s ^ • The precise 

mechanism producing the cut-off is of no relevance here since the collision 
term depends only logarithmically on the value of the cut-off parameter 
The order of magnitude is 


A 





(A. 5) 


For a more detailed account of the shielding mechanism see e.g. TRUB- 
NIKOV (1965). The values of /mA for various temperatures and elec- 
tron densities are tabulated by SPITZER (1962). 


To evaluate the dominating contribution to the collision term which arises 

from scattering by small angles in the range A* © £ A it sufficies to 

expand the quantity f (■?, Iff t) £„(?, t) appearing in (A. 2) in powers of 

® I P , # 

the velocity transfers *u - Tt and v - ^ up to terms of second order. The 
integral over d 9L- may then be worked out with the result 


C, A) * - fa v h* <?>«>*) 


•ft, . 27 A A «*«£ /to* 


(A. 6) 
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This is the desired approximation for tike collision term in the Boltsmann 
equation. 


TRANSPORT EQUATION FOR BULK VELOCITY 


The bulk velocity "v [x, t) is defined by 
a 

*<(*/*) - f&t-fvfa?,*) 

The Boltzmann equation leads to the following transport equation for 



where the pressure tensor P^x, t) is defined by 


(A. 8) 


<!**« UyjVi, (U- */«,*- '£}'f H (*,*,+) (A.< 

Making use of partial integrations the collision integral J* d^u Ja0 may 

be rewritten in the form 


f A £ — fir / A * 5 tWtyMV < A - 10 > 

To proceed with the evaluation of this integral we now assume that the 
partition functions f^Cac, vT, t) may be approximated by Maxwellian distri- 
butions around the bulk velocity v q (x, t) characterized by ^ temperature 
T o (x. t): 


f- ^,(7/; ( >»u/iTkX, J* A // < A - “) 


In this case the pressure tensor takes the well-known form 


2*'*- i'Vfez: 


and the collision integral may be evaluated in closed form with the 


ft* d l u - 1 - . $■ • t ii L &f- Q J.' S r . Z .l I 

' I 7 - ^ I 'iij, l j 


where v a denotes a mean thermal velocity 
op 




and the function G(x) is given by 


G(x) . 


(a. li) 

result 
(A. 13) 


(A. 14) 


(A. 15) 


in terms of the error function ti (x) . 


